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ABSTRACT 

To reliably detect the cosmic microwave background (CMB) anisotropy is 
of great importance in understanding the birth and evolution of the Universe. 
One of the difficulties in CMB experiments is the domination of measured CMB 
anisotropy maps by the Doppler dipole moment from the motion of the antenna 
relative to the CMB. For each measured temperature the expected dipole compo- 
nent has to be calculated separately and then subtracted from the data. A small 
error in dipole direction, antenna pointing direction, sidelobe pickup contami- 
nation, and/or timing synchronism, can raise significant deviation in the dipole 
cleaned CMB temperature. After a full-sky observational scan, the accumulated 
deviations will be structured with a pattern closely correlated to the observation 
pattern with artificial anisotropics on large scales, including artificial quadrupole, 
octopole etc in the final CMB map. Such scan-induced anisotropics on large scales 
can be predicted by the true dipole moment and observational scan scheme. In- 
deed, the expected scan-induced quadrupole pattern of the WMAP mission is 
perfectly in agreement with the published WMAP quadrupole. With the scan 
strategy of the Planck mission, we predict that scan-induced anisotropics will 
also produce an artificially aligned quadrupole. The scan-induced anisotropy is 
a common problem for all sweep missions and, like the foreground emissions, has 
to be removed from observed maps. Without doing so, CMB maps from COBE, 
WMAP, and Planck as well, are not reliable for studying the CMB anisotropy. 

Subject headings: cosmic microwave radiation - cosmology: observations - meth- 
ods: data analysis 
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1. Introduction 

The cosmological principle states that the universe is homogeneous in the large-scale 
average. However, after the Wilkinson Microwave Anisotropy Probe (WMAP) first-year 
release, the quadrupole (with multipole moment / = 2) and octopole {I = 3) patterns of the 
cosm ic microwave background (CM B) map were reported to be aligned along a particular 



axis ( Ide Oliveira-Costa et al. 20041 ) . Later, with more significant statistics the quadrupole 
plane and the three octopole planes are far more aligned, three of these planes are orthogonal 
to the ecliptic, the normals to these planes are aligned with the direction of the dipole and 
with t he equinoxes, and th e remaining octopole plane is orthogonal to the supergalactic 



plane (jSchwartz et al. 2004j); the next two harni onics also appear to be aligned with the 



quadrupole and octopole (jhand &: Magueijo 2005 ). Such odd large-scale anomalies in CMB 



anisotropy, indicating the universe being arrayed around special axes, were also found by 
more works on the following released three-, five- and seven-year WMAP data. 

For CMB experiments and for studies of the very early universe, it is greatly important 
to find out whether the curious behavior in the low / CMB modes being just a coincidence 
produced by primordial fluctuations in our particular Hubble volume (by cosmic variance). 



or caused by a measurement error (IHinshaw et al. 20071 : ICho 20071 ). Here we point out 



that these "axes of evil" do exist due to the nature of CMB observations: the Doppler 
dipole is around the direction of the solar system's motion to the CMB; the inhomogeneous 
integration-time maps of both WMAP and Planck observations have characteristic patterns 
with the ecliptic plane being most sparsely observed and two poles over-sampled; further- 
more, the dipole axis depends on the motions of the antenna to the solar system, the solar 
system to the Galaxy, and the Galaxy to the local supergalaxy, etc. Comparing with the 
CMB anisotropy, the Doppler dipole signal is very strong, thus its predicted amphtude has to 
be calculated and removed from the raw data of any observation during the scanning process, 
observation by observation. An error in the predicted dipole signal will cause a deviation 
in the dipole-cleaned temperature, which will be accumulated into artificial anisotropics in 
the final map after completing a full-sky survey. These artificial anisotropics might carry 
traces of the dipole and scanning scheme; therefore, the combining effect of the dipole and 
scanning scheme has to be carefully inspected and removed in CMB experiments. On the 
other hand, the scan-induced anisotropy in CMB maps might also be useful in monitoring 
motions of our local universe. 

In this paper, we show in §2 that measurement and calculation errors relevant to the 
Doppler dipole can induce considerable deviations in observed CMB temperatures. We 
further demonstrate in §3 that such dipole relevant deviations combining with observational 
scan pattern will produce artificial anisotropics in resultant CMB maps, particularly at 



- 3 - 



the largest angular scales, and show how can such scan-induced pattern on large scales be 
predicted by the true dipole moment and observational scan scheme. In §4 we show that the 
expected scan-induced quadrupole pattern of the WMAP mission is perfectly in agreement 
with the published WMAP quadrupole. With the scan strategy of the Planck mission, we 
find that there will also exist an artificial aligned quadrupole. Finally, in §5 we demonstrate 
that it's possible to remove the observational scan induced effect, and only after doing so 
can the CMB anisotropy be reliably detected. 



2. Temperature Deviation from Doppler Dipole 

2.1. Dipole moment. 

For a sky pixel p with the unit direction vector n[p), the Doppler signal 

T 

T,(p) = ^v-n(p), (1) 

where Tq = 2.725 K is the monopole temperature, v the joint velocity of antenna relative 
to the CMB, c the speed of light. The Doppler signal on all sky pixels gives a dipole 
moment with a pattern of half- hot, half-cold structure along a special axis, which strongly 
dominates observed full-sky microwave maps with amplitude greater than 3 mK, two orders 
of magnitude stronger than the CMB anisotropy. 



2.2. Dipole signal in measured temperature 

Let pit) denote the sky pixel pointed by the antenna at time t, the dipole component 
in the measured temperature can be calculated by 

TMt)) = ^ v(t) ■ n{p{t)) = ^ v'(t) ■ n'ipit)) , (2) 

where v(t) is the spacecraft velocity relative to the CMB rest frame, and n{p{t)) is the unit 
direction, or line-of-sight (LOS), of the antenna at time t in the CMB rest frame. v'(t) 
and n'{p{t)) are corresponding vectors in the spacecraft coordinate system. To remove the 
dipole signal in CMB map, has to be evaluated for each observation by Eq.[2]and then be 
subtracted from the observed time-ordered data. 

An LOS error will produce a pseudo signal in the calculated Td 

ATMt)) = ^ v(t) ■ AnMt)) = ^ v'(t) ■ An; , (3) 
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where An^(p(t)) and An'^ are the LOS error vectors in a sky coordinate system and in the 
spacecraft frame, respectively. In the spacecraft frame the LOS error An'^ is a constant 
vector. After subtracting the calculated dipole signal T^{p{t)), a temperature deviation 
—AT^{p{t)) will be left in the resulting CMB temperature on the sky pixel p{t). 

Even if the real LOS is exactly determined, other observational errors can also produce 
pseudo dipole signal just like an LOS error. For example, we found a 25.6 ms timing offset 
between the spacecraft attitude and radiometer output timestamps in the raw WMAP time- 
ordered data, which can generate a pseudo signal in calculating the Doppler dipole as if there 
were ~ 7' LOS error (Liu, Xiong & Li 2010). The error in the direction of CMB dipole, or 
in the vector v in Eq.|2l can also produce a pseudo dipole signal like a LOS error. The CMB 
dipole determined by the COBE or WMAP mission has a direction error great than ~ 6' 
( iBennett et al. 2003al ). which cannot be ignored in producing pseudo dipole signals. 



2.3. Sidelobe contamination from dipole. 

A microwave antenna has response to a large solid angle, not merely in the main beam 
along the LOS. The complete sidelobe response can be described by a full sky map of 
normalized gain G' in the spacecraft coordinate system, in which the antenna is static and 
G' does not change with time, with the summation of all gains (including the main beam) 
being equal to (the number of pixels in the map). The sidelobe pickup comes mainly from 
the dipole moment, which is the dominating signal source in the foreground-cleaned sky. For 
an observation to the sky pixel p{t) at time t, the recorded signal is contaminated by the 
dipole induced sidelobe pickup 

Ts{p{t)) = J2G'ip's)TM)/N, (4) 

p's 

where the summation goes over all sidelobe pixels p'^ in the spacecraft frame, and Td{p's) is 
the dipole amplitude of the sky pixel ps- The sidelobe pickup Ts{p{t)) calculated by Eq.H] 
also has to be subtracted from the observed temperature. An error AG' in sidelobe response 
will produce a pseudo signal in the sidelobe pickup calculated by Eq.H] 

AT.(p(t)) = ^ v'(t) ■ n'(p;) = ^ v'(t) • An; (5) 

c ^-^ A c 

p's 

where it's self-evident that An^, the result of the summation, is a constant vector. Thus it 
can be seen as an equivalent LOS error. After subtracting the calculated sidelobe pickup, a 
temperature deviation —ATs{p{t)) will be left in the resulting CMB temperature at the sky 
pixel p(t). 



- 5 - 



2.4. Overall LOS error. 

In stead of Eq. [3] and Eq. [5l we can take An' = An'^ + An'^ as an overall LOS error and 

use 

AT(p(t)) = -— v(t) ■ An(p(t)) = -— V(t) ■ An' (6) 
c c 

to estimate the overall deviation in the observed temperature at pixel p{t) after subtracting 

the dipole signal and sidelobe pickup contamination. 



3. Scan Induced Anisotropy 

In Eq.ini the velocity v(t) is determined mainly by the solar system's movement with 
respect to the CMB, and is hence approximately a constant vector in the Galactic coordinate 
system. The overall LOS error An' is a constant vector in the spacecraft frame. During a 
scan period that vectors An{p{t)) are pointing within the hot-half /cold-half area of the 
dipole moment pattern in the Galactic coordinate system (in other words, the dot-product 
v'(t) ■ An' are positive/negative in the spacecraft frame), temperature deviations AT{p{t)) 
will be negative/positive and temperatures T{p{t)) lower/higher, and the corresponding sky 
pixels p{t) will draw a continual cold/hot trace along the scan trajectory in the sky. Finally, 
when a full-sky survey is completed, all deviations AT{p{t)) will be combined into a map. 
The final deviation AT{p) at a certain pixel p can be calculated by 

AT{p) = J2 AT(p(t))/iV, (7) 

p{t)=p 

with Np being the number of observations to the pixel p. The final deviation map is structured 
with a pattern that resembles the exposure map (see Fig. |2l H]). Our IDL source code for 
calculating AT{p) has been made publicly availabl^. 

For a continuous-sweep mission, the design of scanning strategy is constrained by a 
number of technical factors, e.g. solar irradiation, on-board fuel consumption, contact with 
ground antennas etc., it is very difficult to achieve a homogeneous survey. With inhomoge- 
neous coverage, continuous-scan induced anisotropics should consequentially be inhomoge- 
neous with artificial spherical harmonics of low-order moments. In fact, the WMAP team 
already realized that t he large-scale non- G aussian anisotropy features are very similar to the 



WMAP scan pattern (ISpergel et al. 20061 ) . just as expected by the scan- induced anisotropy 



^http: / /dpc. aire. org.cn/ data/ wmap /09072731/release_vl /planck_sim/ 



Fig. 1. — Number of observations per pixel in Galactic coordinates. Simulated one-year scan 
for the Planck mission. Right panel: From the Ql-band one-year scan scheme of the WMAP 
mission. 

analyzed above, and we found that the pixel temperature is significantly correlated with the 
number of observations (jhi et al. 2009 ) . 



4. The WMAP and Planck Missions 
4.1. Planck 



We produce a scan scheme ip(t ) | with the Planck scanning strategy ( jPupac fc Tauber 20051 ) 



and HEALPix ( iGorski et al. 20051 ) resolution parameter Ngide = 1024. The left panel of 
Fig. [1] shows the exposure map for a one-year full-sky survey. Along with the simulated scan 
scheme, we calculate the temperature deviation series {AT(p(t))} by Eq.|6]with an assumed 
overall LOS error An' of 1' and along the X-, Y-, and Z-axis in the spacecraft coordinate sys- 
tem separateljl^. To visually show the trace formation process along with the Planck survey, 
the top panels of Fig. |5] present the traces {AT(p(t))} for a 10-day scanning, where we can 
see temperatures along a scan trajectory being really structured with alternately appearing 
hot and cold. All temperature deviations after one-year full-sky survey are combined by 
Eq.[7]into deviation maps AT^, ATy, and AT^ shown in the middle panels of Fig.[2l where 
we can see the scan-induced deviations closely resembling the Planck exposure time map. 

The bottom panels of Fig.[2]show the quadrupole components, ATg^x, ^Tg^y, and ATg^z, 
of the three deviation maps shown in the middle panels. With these maps, the artificial 
quadrupole from any overall LOS error An' = {Sx,Sy,6z) can be calculated by ATg = 



^In the spacecraft coordinates, the X-axis is parallel to plane of radiators, the Z-axis is the anti-sun 
direction of the spin axis, and the Y-axia is perpendicular to both. 
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Fig. 2. — Expected scan-induced deviation maps for the Planck mission with 1' overall LOS 
error. From left to right: the overall LOS error is along the X, Y, Z-axis of the spacecraft 
coordinate frame, respectively. Top panels: Temperature deviation traces along simulated 
Planck 10-day scan path. Middle panels: Temperature deviations after a full-sky survey. 
Bottom panels: Quadrupole components of the upper panels. All in Galactic coordinates 
and the units are fik. 

6xATq^x + ^y^Tq^y + 6zATq^z. Fig.|3] shows four possible Planck quadrupole patterns with 
{Sx:,Sy,6z) = (1,1,1), (1,1,-1), (1,-1,1), and (1,-1,-1), respectively. It is expected from Fig. [3] 
that, due to scan-induced anisotropies, there may also exist an artificial quadrupole in Planck. 
It may also be aligned with an axis lying in the ecliptic, (like the case shown in the left panel 
of Fig.[3]), or in the Galaxy's plane (the right panel of Fig.|3]), or more possibly between the 
ecliptic and Galaxy's plane. Note that the exact templates of scan-induced anisotropies can 
be calculated when the real Planck data release is available. Thus the amplitudes shown 
here are only qualitative. The real Planck quadrupole depends on the real overall LOS error 
vector. Our calculations are based on the simulated Planck scanning scheme with assumed 
initial conditions of scan geometry. The real pattern of the scan-induced Planck quadrupole 
could hence be some kind of rotation or sign reversion to the ones shown here. 







Fig. 3. — Possible Planck quadruples 

4.2. WMAP 

We use the WMAP Ql-band one-year scan scheme without any mask (because we don't 
use the temperature data) to produce an expected scan-induced deviation map with the 
resolution parameter Nside = 512, the resultant exposure map is shown in the right panel of 
Fig.[TJ Similar to the Planck mission, we calculate temperature deviation maps AT^, ATy 
and ATz shown in the upper panels of Fig. HI where the bottom panels show their quadrupole 
component^. Different to the Planck mission, the deviation components in Y- and Z-axis 
for the WMAP mission have much smaller amplitudes and the three quadrupole components 
are very similar to each other, thus we can predict the scan-induced quadrupole pattern more 
determinately. The prediction is really in agreement with the released WMAP quadrupole. 

In §[2]and §[3l we illustrated that the LOS error, timing offset and sidelobe uncertainties 
can all cause temperature deviations. We find that an overall LOS error with amplitude of 
about 7' can produce most observed WMAP quadrupole (Liu, Xiong & Li 2010). In previous 
works by Roukema (2010) and Liu, Xiong & Li (2011), two different methods have been 
independently used to detect possible timing offset effects in the temperature map and raw 
data respectively, and fully consistent results that support the existence of such effects have 
been obtained. On the other hand, a sidelobe uncertainty of about 10% can also produce 
a similar artificial quadrupole alone. They can take effect simultaneously, and it's hard to 
exactly distinguish them, thus the best way to remove them together is probably by using a 
template-based approach (see §5.2). 

As for the sidelobe uncertainty issue, a possible way is to improve the sidelobe measure- 
ment accuracy to see if the corresponding temperature deviation can be reduced. However, 
we have found that, although the WMAP team have provided sidelobe response data files 



■^The WMAP mission makes measurements with two antennas and records the difference between the two 
antenna temperatures. In this work, we have used two ways to estimate the temperature deviations relevant 
to the dipole issues in the WMAP mission: one by using the same simulation program for Planck (presented 
in footnote 1), and the other by using the real Ql-band one-year scan scheme (ignoring the mask). These 
two ways give consistent results, and the results shown here in Fig H] and Fig U] are by the real scan scheme. 
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Fig. 4. — Predicted scan-induced deviation maps for the WMAP mission with 1' overall LOS 
error using Ql-band one-year real scan scheme. From left to right: the overall LOS error 
is along the X, Y, Z-axis of the spacecraft coordinate frame, respectively. Upper panels: 
Temperature deviations. Lower panels: Quadrupole components of the upper panels. All in 
Galactic coordinates and the units are yuk. 

in each official release, there is no essential improvement except for some minor posterior 
processing adjustment from WMAPl to WMAP7. We have even seen that the Kl-band 
sidelobe response files are exactly the same for all releases. Since the WMAP spacecraft has 
now stopped working, we probably need to wait for Planck for an improved and independent 
sidelobe measurement, so as to run a possible test on the sidelobe uncertainty issues. 

It has been found that the WMAP CMB temperatures are significantly correlated 
with the corresponding number of observations (ILi et al. 20091 ). The scan induced artifi- 
cial anisotropy is probably the most likely explanation to this phenomenon. With the same 
program used to detect the observation number correlation, we have checked the average 
absolute correlation strength between the simulated scan induced artificial anisotropy and 
the corresponding exposure map, which are 0.58 and 0.53 for WMAP and Planck respec- 
tively, strong enough to be able to leave some detectable scan-pattern hke trace amongst 
real CMB, noise and foreground. Therefore, it is possible to check the scan induced artificial 
anisotropy by the correlation method in both missions. If the Planck image quality is really 
very good, then it will be even possible to see the characteristic pattern of scan induced 
artificial anisotropy directly in its result. 
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5. Can CMB Anisotropy be Reliably Detected? 

5.1. Can the overall LOS error be eliminated a priori? 

The equivalent LOS error described in §2.31 is a common problem in CMB experiments, 
because the Doppler dipole moment dominates all observed full-sky CMB maps. The only 
way to eliminate the equivalent LOS error a priori is to know the antenna An response 
exactly, but this is very difficult: If the uncertainty of sidelobe response is about 100%, 
i.e. AG'{p'g) ~ G'{p'g), then with the data ffie of WMAP sidelobe respons^, the amplitude 
6 of the equivalent LOS error induced by sidelobe Doppler pickup (ignore the main lobe 
uncertainty) can be estimated by 

6 ~ arcsin(|An'|/|n'|) = arcsin(| ^ n'(p;)|/|n'|) , (8) 

p's 

which are ~ 50' — 75' for different bands. From Fig.HJwe see that 1' equivalent LOS error 
(requiring sidelobe error < 2%) is enough to produce ~ 1 /iK temperature deviation. How- 
ever, this is for the sidelobe only. The main lobe response is much stronger than sidelobe: 
If we take the main lobe efficiency as 98%, then the overall main lobe response is about 50 
times stronger than the overall sidelobe, and hence < 0.04% main lobe response error can 
probably produce temperature deviation with similar amplitude to 2% sidelobe uncertainty. 
If we believe that the temperature deviation should be less than 1 /xK for a highly reliable 
CMB estimation, then these two percentage data (2% for sidelobe and 0.04% for main lobe) 
roughly describe the required antenna response accuracy. 

For the sidelobe gain map of the first yea r WMAP observati on, the overall calibra- 
tion uncertainty is estimated as large as ~ 30% jBarnes et al. 2003h l^ Although we believe 



that the Planck mission has a much better performance, a requirement like < 2% side- 
lobe uncertainty or < 0.04% main lobe gain uncertainty might still exceed the equipment 
l imitation. Besides th e sidelobe contamination, ~ 6' error in the Doppler dipole direction 



(IBennett et al. 2003al ) also cannot be ignored. 



Therefore, to get a reliable CMB map, the overall LOS error has to be considered and 
the scan-induced artificial anisotropy has to be removed posterior. Without doing so, CMB 



''The data file of WMAP sidelobe response can be found at 
,http://lambda. gsfc.nasa.gov/product/map/dr4/farsidelobeinfo.cfm , 

^This 30% is actually an average level uncertainty, which means the total sidelobe efficiency is 30% 
uncertain. Apparently, this is only part of the overall uncertainty which is relevant to not only the total 
efficiency, but also response on each pixel p'^. 
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maps from sweep missions, including COBE, WMAP, and Planck as well, all can only be 
preliminary and yield unreliable results. 



5.2. Can scan-induced anisotropy be removed posterior? 

To remove the scan-induced anisotropy, the main difficulty is that the overall LOS error 
is unknown. In other words, the error vector An' = 5x An^. + 5y An^ + 5z An^ has three 
unknown coefficients 5xi Sy, and 6^- Is it possible to remove unfavorable signals with unknown 
intensities? The answer is yes. Several techniques can do it. For example, observed sky 
microwave maps are contaminated by diffusion foreground emission generated by different 
diffusion processes with known physical mechanisms but unknown intensities. The WMAP 
team used template fits with the unknown intensities as fitted parameters to remove these 
foreground emission effectively, and also used the ILC method - using the internal linear 
combination of observed maps at differ ent frequencies with estimated weights to suppress 



foreground and noise as far as possible (IBennett et al. 2003bl ). Similar meth ods can also be 



used to diagnos e the pseudo dipole and minimize the scan-induced anisotropy (ILiu &: Li 2010 



Roukema 2010f ) 



For a CMB mission, we produce template maps AT^, ATy and AT^ of scan-induced 
anisotropy from its scan scheme with an assumed value of 6 and letting 6^ = Sy = 6^ = 6. 
From the observed CMB map T* and with the templates maps, we can derive the cleaned 
map 

T = T* -{cxAT, + CyATy + c,AT,), (9) 

where the coefficients Cx, Cy and are determined by minimizing the variance of T. The result 
of template-based removal is not dependent on the amplitude 6 that is assumed for the overall 
LOS error. A larger amplitude leads to a lower fitted coefficient in template fitting so that the 
final correction to the CMB map is the same. To show the ability of scan-induced anisotropy 
removal by template fits, we simulate an observed CMB temperature map with the " synfast" 
routine in the HEALPix software package (available at http://healpix.jpl.nasa.gov) and a 
ACDM model power spectrum with a quadrupole of 1351 /iK^ shown in the left panels of 
Fig.O To the temperature map we add deviations from the WMAP scanning scheme with an 
artificial quadrupole of 6492 /iK^ to produce a simulated observed map, shown in the upper 
middle panel of Fig.O The lower middle panel shows the quadrupole component of the 
simulated observed map with a measured quadrupole of 5768 yuK^. After template removal, 
we get a cleaned map (the upper right panel of Fig. [5]) and its quadrupole component with 
an amplitude of 1186 /iK^ (the lower right panel of Fig.[5]). We can see from Fig.[5]that the 
template removal approach can effectively remove the scan-induced artificial effect from an 
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Fig. 5. — Removing scan-induced anisotropy by template fits. Left panels: Simulated ACDM 
model CMB map. Middle panels: Simulated observed map (ACDM model map + artificial 
anisotropics with WMAP scan scheme. Right panels: After template removal. Upper panels: 
Temperature map. Lower panels: Quadrupole component. All in Galactic coordinates and 
the units of temperature are mK. 

observed map, and in both amplitude and structure pattern, the quadrupole component in 
the cleaned map is more reliable than that without the template removing, even the artificial 
quadrupole (6492 /xK^) is much stronger than the real signal (1351 /xK^). 

5.3. Can largest-scale anisotropies be predicted? 

For a CMB scan mission, template maps representing the distribution of scan-induced 
anisotropies can be produced with the observation scan scheme. Should the artificial anisotropy 
be removed with the templates before releasing a final CMB map? The answer obviously 
should be yes. Unfortunately, in the meantime, many would prefer to accept the odd "axis 
of evil" rather than a cleaned result after necessary correction. They worry about the cor- 
rected quadrupole being too low because the released WMAP quadrupole is already lower 
than expected from the standard cosmology model. In exploring nature, however, to make 
experimental data reliable should be more important than to make it satisfy a theoretical 
requirement. Furthermore, in fact current inflation theories cannot predict the CMB fluc- 
tuation amplitude at the largest angular scale if the time when the scale left the horizon 
was as early as before the reheating period while the primordial density perturbations had 
not occurred. If it is true, it is just beyond, not against the standard model. Therefore, 



- 13 - 



there is no a priori reason to reject a low fluctuation power at the largest scales for the 
sake of the standard model. For an observed full-sky CMB map, scan-induced artificial 
anisotropics should be minimized with predicted templates. More reliable CMB large-scale 
anisotropics cleaned from scan-induced deviations might be a new tool to study the very 
early universe, including the circumstances of the early epoch of inflation before the vacuum 
phase transition, and the start and evolution of the primordial density fluctuations. 
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